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ABSTRACT. In the framework of our investigation to explain the micleosynthetic origin of the 
correlated Ca-Ti-Cr isotopic anomalies in the Ca-Al-rich "FUN" inclusion EK-1-4-1 of the Allende 
meteorite, the nuclear-physics basis in the neutron-rich N = 28 region has been updated by including 
recent experimental data on /3-decay properties and microscopic predictions of neutron-capture cross 
sections. Chargcd-particlc and subsequent r-process calculations within an entropy-based approach were 
performed using a complete reaction network. It is shown that there exist two astrophysical scenarios 
within which the observed isotopic anomalies can be reproduced simultaneously; one at low entropies 
(S ~ 10) which confirms the earlier suggested SN la mechanism, and another at high entropies (5 ~ 150) 
which would be compatible with the neutrino- wind scenario of a SN II. 



P , 1. Introduction 

About two decades ago, G.J. Wasserburg and his group at Caltech published two im- 
portant papers (Lee 1978, Niederer 1980) on the identification of correlated isotopic 
anomalies for Ca, Ti and Cr in peculiar inclusions of the Allende meteorite. These 
highly unusual isotopic compositions were attributed to fractionation (F) and unknown 
nuclear (UN) effects (Lee 1978, Niederer 1980) and therefore designated "FUN" anoma- 
lies. It was immediately concluded that those data might provide a clue to the late-stage 
nuclcosynthetic processes which preceeded formation of the solar nebula. However, all 
astrophysical models existing at that time (standard r-process, generic mini r-process, 
explosive carbon burning, local proton irradiation of solar-system material, nuclear- 
equilibrium approaches) encountered severe difficulties when they tried to reproduce 
the observed anomalies of Allende, in particular those of the EK-1-4-1 inclusion; sum- 
marized by the - slightly acquiescent - statement of the authors: "We acknowledge our 
ignorance of the detailed astrophysical processes... " (Niederer 1980 ) which produced the 
FUN anomalies. 

As already suggested by Wasserburg et al. in 1977, a key to a better description of 
the non-standard abundances may be an improved knowledge of the relevant nuclear- 



physics data around doubly-magic Ca. Following this line, in 1982 Sandler, Koonin and 
Fowler proposed, in their n/3-process, specific nuclear-structure properties, i.e. low-lying 
(30 keV) s-wave neutron-capture resonances, in 46 K(n,7) and 49 Ca(n,7), respectively, 
to replace statistical Hauser-Feshbach (HF) cross sections applied in their reaction net- 
work. These resonances were deemed necessary to enhance the HF rates by an order of 
magnitude resulting in a depletion of the above progenitors of 46 Ca and 49 Ti, otherwise 
being considerably overproduced compared to observation. 

2. Experiments at CERN /ISOLDE 

At this stage, our Mainz-Lund-Garching collaboration became interested in the "EK- 
1-4-1 story". As nuclear chemists / physicists interested in astrophysics, our specific 
approach was to systematically investigate the influence of nuclear-structure proper- 
ties of neutron-rich i 6 S to 19K isotopes on the nucleosynthetic production of their Ca 
and Ti /3-decay daughters. From detailed spectroscopic studies of heavy K isotopes at 
CERN/ISOLDE, rather surprising results had been obtained, for example a sudden, 
strong increase of the /3-delayed neutron decay branch (P D ) from about 1% in TV = 27 
48 K to 86% in N = 28 49 K (Carraz 1982). In 1985, we reported the experimental iden- 
tification of the low-lying s-wave resonance for J 71 " = 3/2~ 49 Ca(n,7) (Ziegert 1985b) 
requested by Sandler et al. (1982), from the "inverse reaction" to neutron capture, i.e. 
the decay mode of /3-delayed neutron emission of 472-ms (J 71 " = CP) 50 K. The J v = 1~ 
state was, however, not found at the energy relevant for n/3-process temperatures of 
T ~ 3 — 5 • 10 8 K, but somewhat higher at 155 keV. Furthermore, subsequent mea- 
surements of the partial decay widths of this state (Ziegert 1985a) have put constraints 
on the Breit-Wigner resonance neutron-capture rate, indicating that it will be roughly 
equal to the HF rate for the above n/3-process temperatures, but may well be larger 
by up to a factor 7 for T > 1 • 10 9 , i.e. for typical r-process temperatures. In any case, 
our experimental data did not support the n/3-explanation of the 49 Ti abundance in 
EK-1-4-1 suggested by Sandler, Koonin and Fowler (1982). 

An extension of this work, using experimental data where available together with 
improved TDA and QRPA model predictions, definitely excluded the n/3-process as a 
possible astrophysical scenario (Ziegert 1985a, Hillebrandt et al. 1986, Kratz et al. 1991) 
and favored explosive He-burning in supernovae (SN) with higher neutron exposures of 
about 5-7T0 -5 mol-cm _3 -s as a likely site for the formation of the correlated Ca-Ti-Cr 
isotopic anomalies in EK-1-4-1. As a consequence of the new nucleosynthesis scenario, 
other - more exotic - progenitor isotopes became the "key" nuclides to produce only 
little Ca and much Ca at the same time, in order to cause the abundance ratio of 
48 Ca/ 46 Ca ~ 250 observed in EK-1-4-1. It was shown that in particular the - at that 
time still unknown - /3-decay half- life of the neutron-magic N — 28 "turning-point" 
nucleus 44 S (where the neutron-capture half- life becomes comparable to the /3-decay 
half-life; Ti/ 2 (n) ~ Ti/ 2 (/3)) would be decisive in defining the neutron-capture process 
duration, and consequently is determining the final amount of stable 46 Ca (Thielemann 
et al. 1990). With shell-model predictions for the decay of spherical 44 S between 1.1 s 
(QRPA) and 0.3 s (TDA), the necessary time scale for the neutron irradiation was 
of the order of 100 ms, and the required abundance of the (a,n) neutron source was 



y a n ~ 2 • 10~ 3 mol-cm -3 . It was well recognized that this might be hard to attain in 
He-zones of massive stars during a SN explosion. Already in 1979, Thielemann et al. 
(1979) found that during explosive processing of a 3a layer of standard preshock com- 
position, maximum abundances of 22 Ne (as (a,n) neutron source) of Ino — 10~ 3 and 
resulting free neutrons of Y n ~ 10~ 9 can be realized for the above requested process 
duration. This was also verified in recalculations with updated reaction rates in 1990 
(Thielemann et al. 1990). Due to the strong influence of the neutron poison 25 Mg pro- 
duced by 22 Ne(a,n), rather the alternative 18 and in particular 13 C neutron sources 
were found to provide the best conditions; however, still at the limit to fulfill the above 
constraints. Nevertheless, the explosive He-burning scenario was not excluded imme- 
diately, because at that time considerable uncertainties still existed in the theoretical 
nuclear-data input of far- unstable progenitor isotopes. Unfortunately - or "fortunately" 
when considering the later, very surprising experimental results presented below - we 
never published this status of the "EK- 1-4-1 story" in a refereed journal... 

3. Experiments at GANIL/LISE 

In order to eliminate at least some of the remaining uncertainties in the nuclear-structure 
properties of potential "turning-point" nuclei south of 48 Ca, in a Mainz-Orsay— Dubna- 
GANIL collaboration we have performed two experiments at the doubly-achromatic 
spectrometer LISE of GANIL. As production mechanism, projectile fragmentation of 
a 60 MeV/u 48 Ca beam was used (Kratz et al. 1991, Sorlin et al. 1993, Kratz et al. 
1995, Bohmer 1996). In these experiments, the gross /3-decay properties Ti/ 2 and P n 
of the very neutron-rich isotopes 43 P, 42 > 44 ' 45 S, 44_46 C1, and 47 Ar could be measured 
for the first time (see Table I). Very surprising results were obtained. Compared to the 
model predictions used in our earlier astrophysical calculations, which - quite naturally 
- assumed spherical shapes for these N ~ 28 near-magic nuclides, the experimental 
T1/2 for most of them (except 44 ' 45 C1) were considerably shorter, by about a factor 3 
for 43 P, 42 S and 47 Ar, a factor 4 for 45 C1, up to a factor 10 for 44 < 45 S. Similarly, most 
of the measured P n values were unexpected when assuming sphericity. We therefore 
suggested that our surprising results indicated an erosion of the N = 28 shell-gap below 
48 Ca resulting in an onset of strong quadrupole deformation (Sorlin 1993, Kratz et al. 
1995). Although initially not predicted by any model of whatever sophistication and 
therefore not accepted in the nuclear-structure community, our "unconventional" ideas 
were rather soon confirmed by both theory and COULEX experiments (see, e.g. Werner 
1994, Scheit 1996, Glasmacher 1998). 

4. Ground-State Shapes from a QRPA Parameter Study 

The main motivation for our shell-model parameter study within the quasi-particle 
random-phase approximation (QRPA) with Folded- Yukawa single-particle (SP) wave 
functions (Moller and Randrup 1990) was to determine the nuclear shape - spherical or 
deformed - of the above N ~ 28 isotopes exclusively from a comparison of their experi- 
mental and calculated gross /3-decay properties Ti/ 2 and P n (Bohmer 1996). These quan- 
tities were calculated as functions of quadrupole deformation, where the e2-parameter 



TABLE I 

QRPA-calculation of Ti/2 and P n with optimized ground-state quadrupole deformation of the 
investigated nuclei. In column 2, the deformation parameter £2 is given which leads to the best 
agreement between the calculated 7\/ 2 and P n values (columns 3 and 4) and the experimental 
results (columns 5 and 6). 



Isotope 


Optimized £2 


Theoretical values 


Experimental values 


Ti /2 [ms] 


Pn[%] 


T 1/2 [ms] 


Pn[%] 


43 p 


-0,175 


70 


100,0 


36±4 


100+10 


42 g 


-0, 300 


1300 


0,9 


579±46 


<4 




+0, 300 


1570 


0,1 






44g 


+0,325 


133 


15,3 


123±10 


18±3 


45 g 


-0, 200 


83 


39,8 


72±5 


54±6 




+0, 125 


90 


36,7 






44 C1 


-0,250 


434 


1,8 


438±19 


<8 


45 CI 


-0, 200 


392 


10,4 


412±39 


38±5 


46 CI 


+0, 100 


204 


57,9 


223±37 


60±9 


47 Ar 


-0, 180 


705 


0,1 


613±197 


<1 




+0, 125 


709 


0,1 







was varied in small steps between -0.35 and +0.35, i.e. between strongly oblate and 
strongly prolate. 

The nuclei considered in this study, have 15 to 18 protons and 26 to 30 neutrons. In 
this low-mass region the single-proton and single-neutron levels are well separated by 
1-5 MeV (see, e.g. the Folded- Yukawa SP-diagrams in Moller 1997). For both protons 
and neutrons, pronounced spherical shell-gaps exist for particle numbers 14, 20 and 
28, and deformed gaps develop around 62 ~ +0.2 for particle numbers 16, 26 and 30 
and around 62 ~ —0.25 for particle number 18, replacing the above spherical shell- 
gaps. In any case, more or less independent of the nuclear shape, the low SP-level 
densities together with the /3-decay selection rules yield only a few Gamow- Teller (GT) 
transitions. In fact, in the A ~ 45 mass region considered here, the quantities Ti/ 2 and 
P n are to a large extent determined by one or two strong GT-transitions in the respective 
energy intervals, whereas the effects of forbidden transitions can be neglected. Thus, the 
expected correlation between the gross /3-decay properties and the nuclear shape should 
be unambiguous, at least to distinguish between sphericity and collectivity. 

In the following, let us excmplarily discuss the results and their relation to our 
astrophysical "turning-point" nucleus 44 S. According to the prediction of most global 
models, the neutron-magic nucleus 44 S should be spherical in its ground state. A QRPA 
calculation in the spherical limit leads to a rather long half-life of Ti/ 2 = 1.12s and a 
maximum possible neutron emission probability of P D = 100% (Sorlin 1993, Kratz et al. 
1995). In this case, the two (lowest-energy) ^7/2 — > 7rf 7 / 2 GT-transitions determining 
both gross /3-decay properties lie above the neutron separation energy, at about 5 MeV 
in the daughter nucleus 44 C1. Much better agreement with the experimental values of 
Ti /2 = (123 ± 10) ms and P n = (18 ± 3) % are obtained if we perform the QRPA cal- 



culation with a prolate deformation of e 2 ~ 0.3. For this shape, GT-transitions between 
deformed Nilsson-orbitals of fj/ 2 shell- model origin are shifted down in energy to near 
the ground state of 44 C1, resulting in a shorter half-life of T 1 / 2 =133 ms together with a 
lower /3-delayed-neutron emission probability of P n = 15.3%. From this comparison, one 
clearly can exclude the spherical shape for 44 S. As mentioned above, in the meantime our 
conclusion has been confirmed by intermediate-energy Coulomb excitation experiments 
of B(E2)-values (Scheit 1996, Glasmacher 1998). 

In the same way, the results on the non-spherical ground-state shapes for the other 
N ~ 28 nuclei were obtained. Table I summarizes the deformation parameters for which 
our QRPA calculations yielded the best agreement with the experimental values. In a 
number of cases, e.g. for 42 ' 45 S and 47 Ar, our approach is not able to determine whether 
these isotopes are prolate or oblate deformed, but can again exclude a spherical shape. 
This general distinction between collectivity and sphericity is, however, of major impor- 
tance as input for improved predictions of neutron-capture cross sections (Rauscher et 
al. 1998) entering our astrophysical network calculations discussed below. 

5. Theoretical Neutron-Capture Rates 

Apart from nuclear masses and /3-decay properties, also neutron-capture cross sections 
(<r n ) are important input data for calculations of non-equilibrium r-process-like nucle- 
osynthesis scenarios. For typical astrophysical neutron energies, two reaction mecha- 
nisms have to be considered (Rauscher 1995): (i) Capture reactions where a compound 
nucleus (CN) is created; in this case cr n (CN) values are commonly calculated with the 
statistical Hauser-Feshbach model, (ii) Reactions where a neutron is directly captured 
into bound states, i.e. direct capture, (DC). The DC-mechanism can dominate the 
neutron-capture for nuclei with low level density in the vicinity of the neutron sepa- 
ration energy S n . This is especially the case for light-mass nuclei, for isotopes close to 
magic neutron numbers, and for nuclei far from /^-stability with low S n values. There- 
fore, in the present study, in addition to CN cross sections (Rauscher et al. 1998) also 
the DC contributions of all even-even (ee) nuclei (because of level-density arguments 
one can assume to first order that DC will be more important for capture on ee targets 
than on odd-nucleon ones) between 40 S to 74 Fe were calculated with the code TEDCA, 
as follows. 

The theoretical cross section a th is given by a sum over each final state i (Kim et al. 
1987) 

a^^CfSiO-r . (1) 

i 

In our case the isospin Clebsch-Gordan coefficients d are equal to unity. The spectro- 
scopic factors Si describe the overlap between the antisymmetrized wave functions of 
target+n and the final state. In the case of one-nucleon capture on even-even deformed 
nuclei, the spectroscopic factor for capture into a state i, which has an occupation 
probability vf in the target, can be reduced to (Glendenning 1983) 



Si = 1 - vt 



(2) 



The corresponding probabilities v\ are found by solving the Lipkin-Nogami pairing 
equations (Moller 1995). 

The factors of c in Eq. |l| are essentially determined by the overlap of the scattering 
wave function in the entrance channel, the bound-state wave function and the multipole- 
transition operator. The potentials needed for the calculation of the before-mentioned 
wave functions are obtained by applying the folding procedure. In this approach, the 
nuclear density of the target px is folded with an energy and density dependent effective 
nucleon-nucleon interaction w e ff (Kobos et al. 1984) 

V(E,R) = XV F (E,R)=xJp T (r)w eff (E,p T ,\R-r\)dv , (3) 

with R being the separation of the centers of mass of the two colliding nuclei. The 
interaction w e tf is only weakly energy dependent in the energy range of interest (Ober- 
hummer and Staudt 1991). The density distributions pt were calculated from the folded- 
Yukawa wave functions. 

The only remaining parameter A was determined by employing a parametrization of 
the volume integral I 

I{E) = -^ Jv F {R,E)R 2 dR , (4) 

expressed in units of MeVfm 3 , and with the mass number A of the target nucleus. 
Recently, the averaged volume integral Io was fitted to a function of mass number A, 
charge Z and neutron number N for a set of specially selected nuclei: (Balogh 1994, 
Krausmann ct al. 1996) 

N — Z 

la = 255.13 + 984.85^~ 1/3 + 9.52 ■ 10 6 . (5) 

Thus, the strength factor A can easily be computed for each nucleus by using 

A=f • (6) 

For the bound states, the parameters A are fixed by the requirement of a correct repro- 
duction of the separation energies. 

The calculated CN and DC cross sections as well as the deformation parameters 
used for the calculation of the single-particle levels are shown in Table II. Experimen- 
tally known neutron separation energies in the final nuclei were taken from a recent 
data compilation (Audi and Wapstra 1993); otherwise predictions from the finite-range 
droplet mass model FRDM (Moller 1995) were used. Furthermore, in this table /3-decay 
half-lives (Ti/ 2 (f3)) are compared to theoretical lifetimes for neutron capture (T^^n)). 
Experimental Ti/ 2 ((3) were again taken from recent compilations (Audi and Wapstra 
1993, Pfeiffer et al. 2000), and - where necessary - theoretical values were obtained by 
using the QRPA code of Moller and Randrup (1990). The Ti/ 2 (n) were computed from 
our present DC+CN results with neutron number densities of 3 • 10 19 cm~ 3 for S and 
Ar target isotopes and 6 • 10 20 cm -3 for Ti, Cr and Fe isotopes, respectively. For the Ca 



TABLE II 

Calculated 30-keV (cm.) Maxwellian averaged neutron-capture cross sections < a >^okeV for 
CN and DC. The column labeled "%" gives the fraction of direct capture in the total cross 
section. Also shown are the deformation parameters £2 and the neutron separation energies S n 
of the final nuclei, target+n. The neutron-capture half- lives Ti/ 2 (n) were computed with the 
values from column "DC+CN" and neutron number densities of 3 • 10 19 for S and Ar target 
isototpes, and 6 ■ 10 20 cm~ 3 for Ti, Cr and Fe isotopes, respectively. The values for the Ca 
isotopes were taken from an earlier calculation using the same model (Krausmann et al. 1996). 



Target 


£2 




DC 


CN 


DC+CN 


% 




1/ j£ \!~' J 






[MeV] 


[mb] 


[mb] 


[mb] 




[s] 


[s] 


4Ug 


+0.24 


3.8238* 


0.4246 


0.0851 


0.5097 


83 


0.218 


8.8+2.2* 


42 g 


+0.30 


3.3114* 


0.9466 


0.0202 


0.9668 


98 


0.115 


0.579+0.046* 


44 g 


+0.30 


1.3344 


0.014 


0.0044 


0.0184 


76 


6.040 


0.123+0.010* 


46 Ar 


-0.18 


4.2590* 


0.5295 


0.1203 


0.6498 


81 


0.171 


8.4+0.6* 


48 Ar 


-0.22 


1.7074 


0.0427 


0.0144 


0.0571 


75 


1.946 


0.11 


50 Ar 


-0.28 


1.0804 


0.0016 


0.0030 


0.0046 


35 


24.15 


0.05 


56 Ti 


+0.13 


2.1936* 


0.0147 


0.1305 


0.1452 


10 


0.038 


0.15+0.03* 


58 Ti 


-0.10 


2.4244 


0.0185 


0.0797 


0.0982 


19 


0.057 


0.047+0.010* 


60 Ti 


-0.02 


2.1194 


0.0165 


0.0403 


0.0568 


29 


0.098 


0.054 


62 Ti 


-0.04 


0.5878 


0.0068 


0.0030 


0.0098 


69 


0.569 


0.018 


64 Ti 


+0.06 


0.4094 


0.0013 


0.0002 


0.0015 


87 


3.561 


0.039 


66 Ti 


+0.14 


0.3914 


0.0009 


0.0002 


0.0011 


82 


4.960 


0.013 


62 Cr 


+0.30 


2.9134 


0.0119 


0.3846 


0.3965 


<1 


0.014 


0.187+0.015* 


64 Cr 


+0.05 


1.8614 


0.0170 


0.0353 


0.0523 


33 


0.106 


0.044+0.012* 


66 Cr 


+0.10 


2.2374 


0.0126 


0.0673 


0.0799 


16 


0.071 


0.071 


68 Cr 


+0.16 


1.8274 


0.0129 


0.0268 


0.0397 


32 


0.140 


0.026 


72 Fe 


+0.14 


2.2474 


0.0101 


0.0572 


0.0673 


15 


0.083 


0.089 


74 Fe 


+0.07 


2.0564 


0.0104 


0.0360 


0.0464 


22 


0.120 


0.052 


76 Fe 


+0.06 


0.0584 


0.0050 


3 • 10" 6 


0.0050 


>99 


1.12 


0.045 



experimental /3-decay half-life (Audi and Wapstra 1993, Pfeiffer et al. 2000) 
experimental S n value (Audi and Wapstra 1993) 



nuclides, the values were calculated earlier by Krausmann et al. (1996) with the same 
model. 

From Table II one can see nicely the importance of direct capture when approaching 
the magic neutron number N = 28 (S, Ar), but also the increasing contribution of DC 
to the total cross section when approaching the neutron drip-line. As already mentioned 
above, "turning-points" are reached in our r-process scenario for those isotopes where 
the half-life for capturing the next neutron becomes longer than the /3-decay lifetime. 
Consequently, the nuclear-physics properties of N = 28 44 S, N = 30 48 Ar, N = 36 58 Ti, 
TV = 40 64 Cr and N = 48 74 Fe are of primary importance. 



6. Astrophysical Network Calculations 

In the late phase of massive stars, prior to the occurrence of rapid neutron-capture 
nucleosynthesis, reactions with charged particles such as (p,7), (p,n), (a,n) and (a, 7), 
as well as their inverse reactions, /3 + -decays and electron captures take place at high 
temperatures of T 9 > 3. The process paths as well as the development of the stars can 
be studied in full network calculations nowadays. The reliability of such calculations 
depends, among other things, also on the sophistication and internal consistency of the 
nuclear-physics input of the involved isotopes. Hence, in our nucleosynthesis calcula- 
tions all nuclear input data, such as masses, neutron-separation energies, level-energies, 
their spin and parity, the ground-state shape, neutron-capture cross sections and decay- 
properties were taken into account. If available, experimental data were used, otherwise 
the input data were obtained from large-scale calculations using modern macroscopic- 
microsocpic models. 

As already pointed out in the introduction, the essential goal of our present investi- 
gation was to find a realistic nucleosynthesis process which reproduces simultaneously 
the complete Ca-Ti-Cr isotopic FUN anomalies observed in the EK-1-4-1 inclusion. 
The required stellar parameters, in particular neutron density (p n ), radiation-entropy 
(S oc Tg I p) and time scale of the process (r), are then expected to give tight constraints 
on the possible astrophysical conditions within the chosen scenario. As in Freiburghaus 
et al. (1999), our calculations follow a hot blob of matter with entropy S: (1) Initially, it 
consists of neutrons, protons and some a-particles in a ratio given by nuclear statistical 
equilibrium (NSE) for a specific proton-to-nucleon ratio Y e , and expands adiabatically 
and cools from Tg ~ 9 to freeze-out conditions of Tg ~ 1. (2) Thereby, the nucleons 
and a-particles combine to heavier nuclei (typically Fe group), with some free neutrons 
and a-particles remaining. (3) For high entropies, an a-rich freeze-out from charged- 
particle reactions occurs for declining temperatures, producing nuclei in the mass range 
A ~ 80-100. (4) Finally, below Tg ~ 3 these nuclei can capture the remaining neutrons 
and undergo r-process-like processes. We consider a simple model of an adiabatically ex- 
panding homogeneous mass zone with an expansion timescale relevant to the supernova 
problem. It is generally accepted now that this approximation is justified by hydrody- 
namic and semi-analytical simulations of the evolution of the high-cntropy wind of a 
SN. The calculations were performed for a grid of entropies (S — 3, 10, 20,. . . , 300 
/ce/baryon) and proton-to-nucleon ratios (Y c = 0.42, 0.45, 0.48), the latter "measuring" 
the neutron-richness of the initial composition. 

Because of the chronological order of the so-called a-process and a neutron-capture 
process in the late phase of a massive star, our computer-code was organized in two steps. 
That part of the network which calculates the charged-particle reactions, was developed 
by F.-K. Thielemann (1990) and extends from iH to 46Pd. On the neutron-rich side 
all nuclei up to the neutron-drip line were included. On the proton-rich side, only the 
first three isotopes away from /3-stability were taken into account. After the initial 
calculations of the a-process, subsequent calculations were performed with the r-process 
code from C. Freiburghaus for those cases where one or more free neutrons per seed- 
nucleus (Y n /Y scc( i > 1) were still existing. A detailed description of the physical processes 
and the necessary system of equations is given in (Freiburghaus 1999). These r-process 
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Fig. 1. Calculated isotope abundance ratios of 4S Ca/ 46 Ca, 48 Ca/ 50 Ti and 48 Ca/ 54 Cr as a func- 
tion of entropy. Here, a proton-to-nucleon ratio of Y c = 0.42 has been chosen. The dashed lines 
represent the respective abundance ratios observed in the EK- 1-4-1 inclusion of the Allende 
meteorite. 

calculations are fully dynamical and are not restricted to a (n,7)-(7,n) equilibrium, as 
were former calculations. 

6.1. Isotope Abundances as a Function of Entropy 

In Fig. 1 the calculated isotope abundance ratios for 48 Ca/ 46 Ca, 48 Ca/ 50 Ti and 48 Ca/ 54 Cr 
at Y e = 0.42 are displayed as a function of entropy (3 < S < 300). The dashed lines 
represent the values observed in EK- 1-4-1; i.e. 48 Ca/ 46 Ca ~ 2 50 , 48 Ca/ 50 Ti ~ 3, and 
48 Ca/ 54 Cr ~ 1. The isotopes 48 Ca, 50 Ti and 54 Cr are the most neutron-rich stable ones 
for each element, and all three are significantly overabundant compared to "normal" 
solar matter. Moreover, even their ratios among each other are enriched in 48 Ca by 




Mass number, A Mass number, A 

Fig. 2. Isotopic abundances calculated for a proton-to-nucleon ratio Y e = 0.42 and 

entropies of S = 10 fce/baryon (left part) and S = 150 fce/baryon (right part). 



factors of 3-4 compared to solar. Therefore, over and above a reproduction of the abun- 
dances of the stable isotopes within each Z-chain, the most sensitive test of our entropy- 
based model will be the above abundance ratios. 

As can be seen from Figure 1, there are two entropy regions for which these iso- 
tope abundance ratios in EK-1-4-1 can be reproduced, one at low entropies (S ~ 10 
/ce/baryon) and the other at high entropies (S ~ 150 fce/baryon). In the first case, 
the Ca-Ti-Cr isotopes are synthesized in an a-rich freeze-out where essentially no free 
neutrons are available for subsequent neutron captures. Hence, the neutron-rich stable 
isotope 48 Ca, for example, is formed directly. This is in agreement with the earlier study 
of Meyer, Krishnan and Clayton (Meyer 1996). In the second case, already some amount 
of free neutrons exists that can start a kind of weak r-process after the charged-particle 
freeze-out. Under such conditions, 48 Ca is mainly generated by /3-decay of its progenitor 
48 Ar. And 44 S, indeed, acts as a "turing-point" nucleus, thus avoiding the production 
of large amounts of 46 S the main progenitor of stable isobaric 46 Ca. 

It may be of further interest to compare the whole abundance distributions produced 
under the above two entropy conditions (see Fig. 2). In principal agreement with the re- 
sults obtained by Freiburghaus et al. (1999), for low entropies sizeable abundances from 
the freeze-out of an a-rich Si-QSE (quasi nuclear statistical equilibrium) are obtained 
only in a limited mass range; i.e. within 45 < A < 95. For higher entropies, abundances 
from the freeze-out of a neutron-rich Fe-QSE extend up to heavier masses. They show 
an increasing pattern from A ~ 40 to A ~ 120, with a sudden drop around A ~ 130 
where the neutron-capture path reaches N = 82 neutron-magic isotopes. There are also 
considerable differences in the absolute abundances of the two scenarios, shown in Fig. 2 
in the units of the solar Si abundance (Si = 10 6 ). At low entropies, the even-even Ca- 
Ti-Cr isotopes of interest are produced at a level of about 10~ 4 , whereas for the higher 
entropies they lie roughly three orders of magnitude lower. 



It may also be of interest to consider abundance ratios of e.g. Ca with heavier 
neutron-rich stable isotopes expected - and, indeed, partly observed in CaAl-rich inclu- 
sions - to be overabundant compared to solar, such as 58 Fe, 64 Ni, 66 Zn and 96 Zr. As 
may already be inferred from a fleeting glance of Fig. 2, there are large differences for 
the low- and high-entropy scenarios. If some tiny amount of EK- 1-4-1 material could 
still be found, those different abundance patterns might motivate further experimental 
studies with todays much more sensitive techniques. To give an example, in case of the 
formation of the Ca-Ti-Cr abundances in a Si-QSE (see, S=10) the abundance ratios 
of Y( 9 °- 96 Zr)/Y( 48 Ca) should lie between 5-1CT 4 and 5-1CT 10 , i.e. negligible, whereas 
for an Fe-QSE followed by a weak r-process (S=150) the above ratios would be or- 
ders of magnitude larger (in detail, 56, 1.7, 11, 45, 120 for 9 °. 91 . 92 . 94 . 96 Zr, respectively). 
Under these latter conditions, the calculated Y( 91 ' 92 Zr)/Y( 90 Zr) isotopic ratios would 
be exactly the solar r-abundance (Y rjQ ) ratios (Kappeler et al. 1989, Arlandini et al. 
1999), whereas the 94 > 96 Zr/ 90 Zr ratios already show considerable overabundances of the 
heaviest Zr isotopes compared to Y r;Q by factors of about 8 and 12, respectively. 

In the context of the question, which of our two nucleosynthesis scenarios - the 
low- or the high-entropy approach - has most likely produced the FUN anomalies in 
EK- 1-4-1, it is interesting to remember that in fact the very first measurements were 
made for Ba, Nd and Sm (Lugmair et al. 1978, McCulloch and Wasserburg 1978a, 
1978b). They found excesses of the unshielded heaviest isotopes, which are generally 
believed to be pure r-process products. Within our entropy-based parameter study, only 
the high-entropy scenario is able to produce isotopes beyond A~130 with measurable 
abundances (Y(Z>56)>1CT 10 ). At the example of the Sm isotopic ratios, relative to 
r-only 149 Sm, in Table III we compare the weighted averages of the experimental values 
obtained by Lugmair et al. (1978) and McCulloch and Wasserburg (1978) for EK-1-4- 
1 with those of the solar-system s- and r-process abundances and our predictions for 
different entropies in the range S=150-230. It is immediately evident that the measured 
EK-1-4-1 values are undoubtedly not of s-process origin; however, they are also not 
fully compatible with the entropy- (neutron-density) mix required to explain the Y r!0 
pattern. Our calculations show that an r-process-like scenario with a limited range of 
entropies is able to reproduce the observed values. When starting with S=150, the 
absolute abundance of 149 Sm is of the order of Y( 149 Sm)~2-10 -13 and seems to be still 
overproduced, whereas the two heaviest Sm isotopes are still underproduced relative 
to EK-1-4-1. This picture changes gradually with increasing entropies. The absolute 
abundances of all Sm nuclides increase up to S~210 (e.g., Y( 149 Sm)~4-l(r 10 for S=170, 
3-10~ 7 for S=190, 1-10 -5 for S=210) and then decrease again. In parallel, with increasing 
entropy more and more r-matcrial is shifted to the heavier Sm isotopes. This is clearly 
reflected by the abundance ratios Y( 152 ' 154 Sm)/Y( 149 Sm) in the table. Not only when 
considering the ratios, but also when using the absolute Sm abundances best agreement 
with EK-1-4-1 is obtained around S~190-200. These are conditions just at the upper 
end of the high-entropy region where also the correlated Ca-Ti-Cr anomalies in this 
inclusion could be reproduced. It is worth to be mentioned in this context, that a quite 
similar picture is observed for the Nd isotopic abundance ratios measured by Lugmair 
et al. (1978) and McCulloch and Wasserburg (1978 a). A more detailed exploration of 
the astrophysical conditions to produce those isotopes in EK-1-4-1 and other meteoritic 



TABLE III 

Comparison of Sm isotopic abundance ratios observed in EK-1-4-1 (Lugmair et al. 1978, Mc- 
Culloch and Wasserburg 1978b) with solar-system s- and r-process abundances, and with our 
network calculations for different entropies. 



Isotope 


Y s ,0 


Yr,0 


YeK-1-4-1 




Network calculations 




abund.-ratio 


a) 


b) 


c) 


S=150 


S=170 


S=190 


S=210 


S=230 


i4V Sm/ 14y Sm 


1.92 


1.01 


0.99±0.10 


3.05 


1.82 


1.23 


0.93 


0.87 


152 Sm/ 149 Sm 


3.50 


1.71 


0.74±0.10 


0.25 


0.22 


0.77 


2.18 


2.18 


154 Sm/ 149 Sm 


0.088 


1.87 


l.OOiO.ll 


0.035 


0.16 


0.47 


1.54 


1.81 



a) Kappeler et al. 1990 

b) Arlandini et al. 1999 

c) Lugmair et al. 1978, McCulloch and Wasserburg 1978b 



inclusions is presently underway and will be published elsewhere. 
6.2. Possible Astrophysical Scenarios 

We have shown that resultant abundances of neutron- rich stable isotopes in the 48 Ca to 
66 Zn region can be produced in matter expanding from high temperature and density 
for two different entropy regions. In agreement with Meyer et al. (Meyer 1996), we find 
that for low entropies (S ~ 10 fce/baryon) these isotopes are synthesized in the a- 
rich freeze-out from a Si-QSE. In this case, the key nucleus 48 Ca is synthesized directly 
without significant neutron captures. Such conditions are found in Type la supernovae 
(SN la), which are certainly responsible for the production of the "bulk" 48 Ca material 
in our solar system, as outlined by Meyer, Krishnan and Clayton. 

We have furthermore shown that the Ca-Ti-Cr anomalies in EK-1-4-1 can also be 
reproduced with higher entropies of S~150-190 /ce/baryon, however with considerably 
lower yields than with low entropies. Under these conditions, neutron-rich matter is 
formed by the freeze-out of an a-rich Fc-QSE with subsequent r-process-like neutron 
captures. In this case, the abundance distribution is much larger and ranges up to the 
rare-earth region. A realistic astrophysical site for such a scenario would be the high- 
entropy bubble at the outer core of a type II supernova (SN II). Meyer et al. (1996) 
conclude from their detailed NSE network calculations that 48 Ca does not survive such 
an environment, because the material freezes out with too few nuclei at low masses. In 
their model, the QSE nearly exclusively produces nuclei heavier than 48 Ca; therefore, 
they make the rather stong statement: "We therefore rule out this Type II site as an 
origin for 48 Ca, even in matter having the neutron richness of 48 Ca.". Here, we disagree 
with Meyer et al. (1996). In our network calculations, maybe due to the considerably 
improved nuclear-physics data input for neutron-rich isotopes, we still obtain measurable 
quantities of nuclides in the Ca region simultaneously with overabundances of the 
heaviest stable isotopes up to Nd and Sm, also observed in EK-1-4-1 (Lugmair et al. 
1978, McCulloch and Wasserburg 1978a, 1978b). Typical yields for low-mass isotopes 



with S=150 fc B /baryon are 1 • l(r 10 M Q for 46 Ca, 4 • l(n 8 M Q for 48 Ca, 2 • 1(T 9 M for 
50 Ti and 3 • 1CT 8 M for 54 Cr. This is about 2 to 3 orders of magnitude less than in 
SN la. Moreover, SN II events occur approximately ten times less than SN la. Taken 
together, the contribution of synthesized and ejected matter from SN II is a factor 
10 3 to 10 4 less than the "bulk" material from SN la. Hence, we absolutely agree with 
Meyer et al. (1996), that SN II do not produce enough low-mass material to generate 
the solar amounts of 48 Ca, 50 Ti and 54 Cr. However, the amount of matter ejected by 
Type II supernovae may well be sufficient to explain the observations in the specific 
SN condensate EK-1-4-1. In an attempt to explain the correlation of overabundances of 
neutron-rich isotopes over a wide range of elements in CaAl-rich inclusions, Meyer et al. 
(1996) argue that this material, which was formed in the solar accretion disk, contains 
"components not only from many Type la supernovae but from Type II supernovae as 
well.". In our high-entropy approach, however, all FUN anomalies, from 48 Ca up to 
154 Sm observed in EK-1-4-1, arc produced simultaneously in the high-entropy bubble 
scenario of a SN II. Therefore, extending the early notation of Wasserburg (Lee 1978, 
Nicdcrcr 1980), we would now speak not only about correlated Ca-Ti anomalies, but 
about correlated overabundances of neutron-rich Ca-Ti-Cr-Fe-Ni-Zn. . .Zr. . .Ba-Nd- 
Sm isotopes in EK-1-4-1. 

In full agreement with Meyer, Krishnan and Clayton we conclude that certainly 
much work remains to be done, although we have also learned a great deal. It will 
therefore be interesting to see how our insights into FUN syntheses will fit in with 
future astrophysical models. 
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